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If light sterile neutrinos exist and mix with the active neutrino flavors, this mixing will 
affect the propagation of high-energy neutrinos from dark matter annihilation in the Sun. 
In particular, new Mikheyev-Smirnov-Wolfenstein resonances can occur, leading to almost 
complete conversion of some active neutrino flavors into sterile states. We demonstrate how 
this can weaken IceCube limits on neutrino capture and annihilation in the Sun and how 
potential future conflicts between IceCube constraints and direct detection or collider data 
might be resolved by invoking sterile neutrinos. We also point out that, if the dark matter- 
nucleon scattering cross section and the allowed annihilation channels are precisely measured 
in direct detection and collider experiments in the future, IceCube can be used to constrain 
sterile neutrino models using neutrinos from the dark matter annihilation. 



I. INTRODUCTION 



The hunt for dark matter is currently at a very exciting, but also somewhat confusing stage. 
Many unexpected experimental results that have been reported over the past few years can be 
interpreted in terms of dark matter, but all of them could have more mundane explanations, and 
moreover the dark matter interpretations of different experimental data sets do not fit together in 
many cases. For instance, the signals reported by CoGeNT [ll|2], DAMA [3| |4| and CRESST [1] 
appear to be in some tension with the null results from other direct detection experiment, in 
particular CDMS [6( [7] and XENON-100 [8j (see, however, [9l[l0]). Moreover, the dark matter 
parameter regions favored by CoGeNT, DAMA and CRESST do not coincide under standard 
assumptions on the dark matter halo [TTHTG] (see, however, refs. [17-20J). Also, if the recently 
observed anomalies in the cosmic electron and positron spectra [.214 22\) are due to dark matter 
annihilation or decay, this would imply dark matter masses of order 1 TeV (see, for instance, 
[23]). whereas the CoGeNT, DAMA and CRESST hints would indicate dark matter masses of 
order 10 GeV. It is thus clear that many of the potential hints for dark matter must have other 
explanations, and this illustrates that a single experiment might never be able to unambiguously 
identify dark matter. Only matching detections by several different experiments would convince 
the community at large that dark matter has been observed. Fortunately, the toolbox for dark 
matter search is quite large: Direct detection experiments like CoGeNT, DAMA, CRESST, CDMS 
and XENON-100 search for dark matter recoils on atomic nuclei; collider searches at the Tevatron 
and the LHC aim to directly produce dark matter particles and detect them through missing 
energy signatures; indirect searches look for the annihilation or decay products of astrophysical 
dark matter. Among the possible messengers are electrons and positrons, anti-protons, gamma 
rays, and neutrinos. 

A special role is played by searches for neutrinos from dark matter annihilation in the Sun, 
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which are carried out by the Super-Kamiokande [SUES] and IceCube [26l|22] collaborations. Even 
though these searches probe the products of dark matter annihilation, the expected event rates 
are usually determined by the dark matter capture rate in the Sun and thus by the dark matter- 
nucleus scattering cross section. Therefore, these searches, even though indirect, are sensitive to 
the same observables as direct detection experiments and can directly test any potential direct 
detection signal (provided that dark matter can annihilate and that its annihilation products 
include neutrinos). In particular, many astrophysical uncertainties, for instance those associated 
with the local dark matter density, affect the Super-Kamiokande and IceCube searches in the same 
way as the direct searches, making the comparison between those experiments quite robust with 
respect to astrophysics. 

On the other hand, neutrinos from dark matter annihilation in the Sun are strongly affected 
by neutrino oscillation physics. In this paper, we will investigate how the oscillations pattern of 
high-energy neutrinos from dark matter annihilation in the Sun can be modified by the existence 
of sterile neutrinos. Our study is motivated by the results of the LSND [28j and MiniBooNE [29] 
experiments, as well as the reactor antineutrino anomaly [30H32] , all of which can be interpreted 
as hints for the existence of sterile neutrinos with masses of order 1 eV [331435] . (Note, however, 
that even models with two sterile neutrinos cannot resolve all tension in the global data set.) We 
will argue that, if sterile neutrinos exist, neutrinos from dark matter annihilation can encounter 
new Mikheyev-Smirnov-Wolfenstein (MSW) resonances when propagating out of the Sun, and that 
these resonances can potentially convert a large fraction of them into undetectable sterile states. 
This can weaken constraints on dark matter annihilation in the Sun significantly. (The existence 
of new MSW resonances in the presence of sterile neutrinos has also been investigated recently 
in the context of IceCube atmospheric neutrino data ^361 137].) On the positive side, if the dark 
matter-nucleon scattering cross section and the dark matter annihilation channels are precisely 
determined elsewhere, for instance in direct detection and collider experiments, IceCube can be 
used as a sensitive tool for constraining sterile neutrino models. 

The outline of the paper is as follows: In section [TT| we review the relevant aspects of the formal- 
ism of neutrino oscillations and discuss the effect of MSW resonances on the oscillation probabilities 
of high-energy neutrinos in the Sun. We then describe in section III how we compute the expected 
neutrino signal from dark matter annihilation in the IceCube detector, and in section [IV] we show 
how the existence of sterile neutrinos modifies the dark matter constraints from IceCube. We will 
discuss our results and conclude in section IVII 



II. NEUTRINO OSCILLATIONS AND NEUTRINO INTERACTIONS IN THE SUN 

Neutrinos from dark matter annihilation in the Sun probe a very unique regime of neutrino 
oscillations: They are produced in a region of very high matter density (^ 150 g/cm^) at the 
center of the Sun [^, but with energies that can be much higher than those at which neutrino 
oscillations in the Sun are usually studied. 

In the standard three-flavor oscillation framework, it is well known from the study of low-energy 
((9(MeV)) solar neutrinos that strong transitions between electron neutrinos, z/g, and muon/tau 
neutrinos, z/^, z/^, take place in a region where the number density of electrons reaches a critical 
value, given by the Mikheyev-Smirnov-Wolfenstein (MSW) resonance condition (391442] 



Here, E is the neutrino energy, Gi? is the Fermi constant, 9i2 and Am2i are the usual solar neutrino 
mixing parameters, and acp = 1 (— 1) for neutrinos (antineutrinos). The MSW resonance condition 
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can be understood if we recall that according to the Fermi theory of weak interactions the local 
matter potential due to W exchange with an electron, which is felt by electron-neutrinos but not 
by muon and tau neutrinos, is given by \/2acpGi?ne(r), with ne(r) = (e^^e) the electron number 
density at a distance r from the center of the Sun. At high matter density near the center, the 
flavor- diagonal MSW potential is larger than the flavor-off-diagonal neutrino mass term /S.m\-^/2E 
for multi-MeV neutrinos. Thus mixing between and z/^, Vr is suppressed, and mass and flavor 
eigenstate almost coincide. For instance, the flavor eigenstate z/g is almost equal to the mass 
eigenstate V2 for multi-MeV neutrinos produced at the center of the Sun. At low matter density in 
the outer layers of the Sun, on the other hand, the mass terms dominate over the potential term, 
so that the effective mixing matrix is close to the vacuum mixing matrix, according to which 
is mostly composed of z/i. If the change in the matter density is not too fast, neutrinos cannot 
"jump" from one mass eigenstate to another, so that a neutrino produced as an almost pure V2 
will still be in an almost pure V2 state when it exits the Sun. However, its flavor composition has 
changed dramatically, and in fact, the z/g admixture to 1^2 in vacuum is given |C/e2p — sin^ ^12 — 0.31 
(using the standard parameterization [42j and the current best flt values |l3l for the leptonic 
mixing matrix). Thus, almost 70% of the neutrinos are converted to z/^, z/^- on their way out of the 
Sun. The flavor-conversion happens predominantly at the transition between the matter potential- 
dominated and the mass mixing-dominated regime, where the two terms are of similar magnitude. 
This requirement leads precisely to the condition ([T]). 

For energies above ^ 100 MeV (not accessible with conventional solar neutrinos), a second 
MSW resonance appears at a higher density 

A'«- = <.CPC0.«,3^-^. (2) 

This second resonance leads to strong ^ ly/j^, transitions if Amg^^ > 0, and to strong ^ z>^, 
transitions for Amg^^ < 0. 

The requirement that the change in matter density be not too fast (see above) can be made 
more precise. One can show that resonant flavor transitions in the (zj)-sector cease when the 
adiabaticity condition [42j 

fAml \2 1 
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is no longer fulfllled. Here, 7^ is called the adiabaticity parameter at the resonance and \V 
denotes the gradient of the MSW potential V = \/2Gi?A^e at the location of the resonance. Loss 
of adiabaticity thus occurs for small mixing angles, small Am? and high energies. In the case of 
the resonance in the (12)-sector, which turned out to be the solution to the long-standing solar 
neutrino problem, we expect adiabatic transitions below ^ 10 GeV, and non-adiabatic behavior 
above. (Note that, if flavor transitions of solar neutrinos were non-adiabatic, an initial z/g would 
leave the Sun not as a z/2 mass eigenstate, but as a superposition of the form U^i\ui) + U^2^'^^W'2)^ 
with the oscillation phase 0. After averaging over 0, this would lead to a z/g survival probability 
given by 1 — ^ sin^ 2^12, in conflict with the experimental data on solar neutrinos.) 

The neutrino oscillation probabilities in the Sun in the standard three-flavor framework are 
plotted as a function of energy in flgures[l]and[2](black lines). The transition between the adiabatic 
and non-adiabatic regimes at energies around 10 GeV is clearly visible. At typical solar neutrino 
energies of few MeV, the z/g survival probability has the expected value of sin^ 612 — 0.3, while in 
the non-adiabatic regime, it is 1 — ^ sin^ 20i2 — 0.6. (Small deviations from these values can arise 
from the inclusion of three-flavor effects, in particular a non-zero ^13.) 
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Figure 1: Flavor transition probabilities in the Sun as a function of energy for an initial z^e (left), an initial 
(center), and an initial Vr (right). The top plots are for neutrinos, the ones at the bottom are for anti- 
neutrinos. Black lines are for standard three-flavor oscillation, whereas red lines are for a representative 
"3 + 2" model with two sterile neutrinos (see text for details). Absorption and r regeneration effects are 
neglected in these plots. Note that the black dotted lines {vx in the SM) and the black dot-dashed 

lines ^ ^/i in the SM) lie on top of each other since y^,-yr mixing is assumed to be maximal. 



If sterile neutrinos exist, the oscillation phenomenology becomes much richer. Even if vacuum 
oscillations between active and sterile neutrino flavors are negligible because of small mixing angles, 
active-sterile oscillations in matter can be significant, in particular at the high energies relevant 
to neutrinos from dark matter annihilation. The n-flavor MSW potential has the form 

V^{n,- nn/2, -nn/2, -n,/2, 0, . . . ) , (4) 

where the terms containing the neutron density Un originate from coherent forward scattering 
through exchange. These terms are usually neglected in the three-flavor framework since they 
are flavor-universal and therefore cannot contribute to oscillations among active neutrinos. How- 
ever, they become relevant in the presence of sterile states. In particular, there will be additional 
MSW resonances whenever any of the matter potential terms becomes equal to any of the mass 
terms in the Hamiltonian. These MSW resonance can lead to nearly complete conversion of certain 
neutrino (or antineutrino) flavors into sterile states on the way out of the Sun. 

To illustrate this observation, which is the main topic of this paper, we consider a sterile neutrino 
scenario similar to the one that has been shown in Ref. [33j to provide a reasonably good fit to the 
global neutrino data, including the anomalous LSND and MiniBooNE results. The model has two 
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sterile neutrino flavors Ugi, Vs2 and two new mass eigenstates 1^4, 1/5 with mixing parameters 

sin^ 6x2 = 0.32 sin^ 26I13 = 3 x 10"^ sin^ 62^ = 0.45 
Amii = 7.6 X lO"'^ eV^ Am^i = 2.38 x 10"^ eV^ 
Am^i = 0.90 eV^ Amii = 0.47 eV^ 
sin^ 26*14 = 0.086 sin^ 26*15 = 0.060 sin^ 26*24 = 0.088 
sin^ 26*25 = 0.055 sin^ 26*34 = 0.002 sin^ 26*35 = 0.000 
(5i3 = 1.477r ^14 = 1.0867r ^15 = 0.777r 

Here, we use the parameterization 



(5) 



(6) 



for the leptonic mixing matrix, where Rij denotes a rotation matrix in the {ij) plane with rotation 
angle 9ij^ and Rfj denotes a similar rotation matrix which in addition carries a complex phase Sij: 
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In a "3 + 2" scenario like equation ([5]), the new MSW resonances converting active neutrinos into 
sterile ones affect antineutrinos more strongly than neutrinos, but since neutrino cross sections are 
larger than antineutrino cross sections, we expect the impact of sterile neutrinos on dark matter 
searches to be only moderate, especially in detectors like IceCube and Super-Kamiokande which 
cannot distinguish neutrinos from antineutrinos. (Below, we will also discuss a 3 + 3 toy model in 
which effects are larger.) 

The neutrino oscillation probabilities in the Sun for this sterile neutrino scenario are shown in 
figures [l] as red curves. The most striking feature is the strong conversion of (and to some degree 
also Ur) into sterile neutrinos at energies above ^ 200 GeV. Indeed, we can see from equation [l] 
(with the replacements Ou Ou ^ 0, Amli Amli 2^-1 eV^, and -Nn/2) that 

above E ^ 100 GeV, the MSW resonance between active and sterile neutrinos lies within the 
Sun. Therefore, high energy t?^ and i>r produced from dark matter annihilation at the center 
of the Sun will be almost fully converted into sterile neutrinos, leaving as detectable states only 
neutrinos, and antineutrinos from the z/g component of the primary flux. For a given dark matter 
mass, annihilation channel and annihilation cross section, the expected event number in a neutrino 
detector is thus reduced, so that experimental constraints on dark matter annihilation in the Sun 
become weaker. 

In addition to the 3 + 2 scenario, we are also going to consider a 3 + 3 toy model with 3 sterile 
neutrinos. The oscillation parameters in this model are chosen such that each active neutrino flavor 
eigenstate mixes with only one of the sterile neutrinos. This can be achieved by choosing mass 
squared difference Ami = Amli 2^ ^^52 — ^^63 mixing angles 6s = 6^14 — O25 — 0^6 (all 
other active-sterile mixing angles are zero), so that the sterile neutrino sector is a mirror image of 
the active neutrino sector as far as vacuum oscillations are concerned. (Similar models have been 
considered in [^5].) The parameterization of the leptonic mixing matrix is here 



^^3+3 = R36R25RUR2SR13R12 



(8) 
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Unless specified otherwise, we choose Amli — 0.1 eV^ and sin^ 2^^ — 0.03. In general, if Amli, 
Am|2 , AttIq^ ^ Am2i , | Am^i \ , conversions of active neutrinos into sterile neutrinos can be under- 
stood in a simple two-flavor framework as long as the distance travelled by the neutrinos is much 
shorter than the active neutrino oscillation lengths ^21*^ — ^T^E/Am^i and Lf^ — AttE /\/S.rn^^\. 
This remains true even in matter. In this case, the effective two-flavor oscillations between an 
active flavor and its corresponding sterile flavor are affected by an MSW resonance. The resonance 
conditions are, in analogy to equations ([l]) and Q: 



Ne = OOP COS Bu 
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2E v^Gf 
Ami2 1 
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(11) 



We see from these equations that the resonance between z/g and the first sterile flavor eigenstate 
Z/5I will be in the neutrino sector (acp = 1), whereas the z/^ ^ and ^ i^sS resonances 
affect antineutrinos (acp = — 1). This behavior is reflected in figure [2j where we show the flavor 
transition probabilities for all oscillation channels in the 3 + 3 model as a function of energy. We 
see that in a large energy range zvg, i^/j, and are almost fully converted into sterile states. We 
expect that this will lead to a considerable weakening of the limits IceCube can set on dark matter 
capture and annihilation in the Sun. 

Note that this weakening could be even more pronounced if the mostly active neutrino mass 
eigenstates were heavier than the mostly sterile ones, since in that case the MSW resonances for 
second and third generation neutrinos would move from the antineutrino sector to the neutrino 
sector, which is more important for IceCube's dark matter search because neutrino interaction cross 
sections are about a factor of 3 larger than antineutrino cross sections. We do not consider this 
possibility here since relatively heavy active neutrinos would be in potential conflict with cosmol- 
ogy |l6l - [i9] . (These conflict can potentially be avoided in non-minimal cosmologies [481139] and in 
models where the relic abundance of sterile neutrinos is reduced, see for instance references [501 IST] 
for a discussion of such models.) 

Apart from oscillation, the propagation of high-energy neutrinos through the Sun is also affected 
by non-coherent neutral current (CC) and charged current (CC) interactions. NC interactions 
change the neutrino energy, whereas CC interactions lead to absorption and possible reemission of 
neutrinos in the decay of secondary /i or r leptons. Since secondary muons are usually thermalized 
before they decay, reemission of high- energy neutrinos is only possible in the case ofur+X r-\-X^ 
CC interactions ("r regeneration"). In figure |3] we plot the non-interaction ("survival") probability 
for neutrinos from dark matter annihilation on their way out of the Sun as a function of the neutrino 
energy. 



III. SIMULATION TECHNIQUES 

To estimate quantitatively how existing limits on dark matter annihilation in the Sun are 
modified in the presence of sterile neutrinos, we have carried out numerical simulations. We 
compute the dark matter capture rate as a function of the dark matter mass and scattering cross 
section using the formulae from [52j and assuming a local WIMP density of 0.3 GeV/cm^ with an 
isothermal velocity distribution and velocity dispersion 220 km/sec. We assume the annihilation 
cross section to be large enough for the capture and annihilation reactions to be in equilibrium. 
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Figure 2: Flavor transition probabilities in the Sun as a function of energy for an initial z^e (left), an initial 
(center), and an initial Vr (right). The top plots are for neutrinos, the ones at the bottom are for anti- 
neutrinos. Black lines are for standard three-flavor oscillation, whereas red lines are for a "3 + 3" toy model 
with three sterile neutrinos (see text for details). Absorption and r regeneration effects are neglected in 
these plots. Note that the black dotted lines {i^x in the SM) and the black dot-dashed lines 

in the SM) lie on top of each other since v^,-yr mixing is assumed to be maximal. 




^^[GeV] ^^[GeV] ^^^[GeV] 



Figure 3: Survival (= non-interaction) probabilities for high-energy neutrinos from dark matter annihilation 
on their way out of the Sun. We show results for standard three-flavor oscillations, for the best fitting 3 + 2 
model, and for a 3 + 3 toy model. The features in the 3 + 2 and 3 + 3 curves are due to the interplay of 
active-sterile conversion and active neutrino interactions. 



so that the annihilation rate is equal to half the capture rate. We use initial neutrino spectra 
from [53j, which were generated using WimpSim [54j. 

To propagate the neutrinos out of the Sun, we use our own Monte Carlo code, which is capable 
of working with an arbitrary number of neutrino flavors n, and simulates n-flavor oscillations in 
matter as well as NC and CC neutrino scattering in the Sun, including r regeneration. We use the 
nusigma package [531 [5l] to calculate the neutrino cross sections, and TAUOLA [55j for decaying 
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secondary r's. 

In practice, we proceed as follows: We propagate the n-component neutrino state vector ^(t) out 
of the Sun using the rkf 45 Runge-Kutta-Fehlberg algorithm from the GNU Scientific Library [56j 
to solve the evolution equation 



/o 



21 



f/V(t) + V2Gf 



Nn(t) 

2 



N„{t) 
2 





(12) 



After each Runge-Kutta step, we determine randomly if the neutrino undergoes an incoherent 
interactions during that step. The probability for a CC or NC interaction is given by Pcc/NC — 
^CC/Nc/(^CC + ctnc) X [l ~ ( ~ ^(^) (^CC + ctnc))] i where n(r) is the local nucleon number 
density, crQc(NC) is the charged current (neutral current) neutrino-nucleon scattering cross section, 
and Ar is the current Runge-Kutta step size. If it is determined that the neutrino interacts through 
a neutral current, its energy after the interaction is picked randomly from the final state energy 
spectrum calculated using nusigma [53l[54]. Since we are treating neutrino propagation as a one- 
dimensional problem, we assume that the direction of travel does not change, and we continue to 
propagate the neutrino radially outward. In the case of a or z/^ charged current interactions, 
we simply discard the neutrino. In a charged current interactions, the original neutrino is also 
absorbed, but since the secondary r lepton (unlike a secondary muon from a z/^ interaction) decays 
before it is stopped in matter, new high-energy neutrinos can be produced from its decay ("r 
regeneration"). We use TAUOLA [55] to simulate r decay, and propagate all secondary high-energy 
neutrinos out of the Sun individually. 

We compute the expected event rate in the IceCube detector by multiplying the differential 
muon neutrino and antineutrino fluxes at the Earth the effective detector area AQfi{E) [57j and 
then integrating over energy. We have checked that Earth shadowing effects [58l [59] are negligible 
for our results. Note that the effective area given in [57j has been computed from a simulation of 
the full 86-string IceCube detector, whereas the latest published dark matter limits from IceCube 
are based on data taken in the 40-string IceCube configuration and in the older AMANDA-II 
detector. Since we will ultimately use our simulation only to compute ratios of event rates between 
different oscillation models, we expect the systematic bias introduced that way to be small. Note 
also that AQfi{E) as given in |57j is the combined effective area for neutrinos and antineutrinos. 
Since in sterile neutrino models, the relative importance of neutrinos and antineutrinos in the 
IceCube signal changes, we need seperate effective areas for neutrinos {A^^{E)) and antineutrinos 
{A^^{E)). We obtain them according to 



(13) 
(14) 



with the charged current neutrino-nucleon (antineutrino-nucleon) cross section aQQ{E) (ctqq (£")), 
and the muon (antimuon) range d^- (E"^) ((i^+ (£^^)). For simplicity, we assume a one-to-one relation 
between the neutrino energy E and the secondary muon energy E/^: E"^ = (1 — ycc{Ey))Ey^ where 
ycc is the mean charged current inelasticity parameter [60j. We have checked that using full 
differential cross sections would not significantly change our results. 
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Figure 4: Predicted neutrino fluxes (left) and antineutrino fluxes (right) for annihilation of a 1 TeV WIMP 
into in the Sun. We show the total neutrino flux at production, as well as the muon neutrino flux 

at the Earth. For illustration, we also show the flux of secondary neutrinos from r regeneration, as well as 
the flux obtained using the simplified calculation that neglects regeneration and partial energy loss (see text 
for details). Results for standard oscillations are shown in black, results for the 3 + 3 toy model introduced 
in section ini are shown in red. 

We have verified our Monte Carlo code by comparing its predictions to published results 
from [531 [511 EH [62]. 

While the advantage of the Monte Carlo technique is certainly its flexibility, it is also quite 
computationally intensive. Since we are also interested in carrying out parameter scans over differ- 
ent sets of sterile neutrino parameters (see section [v| below) , we have also developed a faster code, 
which does not take into account r regeneration and energy loss in neutral current interactions. 
Instead, it simply considers all neutrinos that interact in the Sun in any way (NC or CC) to be lost 
to detection. Thus, for each given set of oscillation parameters and for each neutrino energy, we 
need to solve the equation of motion only once to determine the oscillation probabilities for those 
neutrinos which do not interact. At each Runge-Kutta step, we also keep track of the interaction 
probability to obtain simultaneously the fraction of neutrinos at the considered energy which leave 
the Sun without interacting. 

We compare the results of our full Monte Carlo simulation to those of the simplified method 
in figure [Ij We also show the flux of secondary neutrinos from r regeneration, and we notice that 
these neutrinos account for most of the difference between the MC results and the ones from the 
simplified method. (Another small contribution to this difference comes from neutrinos that have 
undergone NC scattering, but are still within the accessible energy range.) This conclusion is the 
same for standard three-flavor oscillations (black and gray curves in figure [4]) and for the 3 + 3 
model (red curves). 
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IV. MODIFIED ICECUBE LIMITS ON DARK MATTER CAPTURE IN THE SUN 

In figure [5] we show how the IceCube hmits on spin- dependent dark matter-proton scattering 
need to be modified if sterile neutrinos exist (black and gray lines). For comparison we also show as 
colored lines limits from a number of direct dark matter searches. Solid black lines in figure [5] are 
the pubhshed IceCube hmits from [271 [57]; Dashed and dotted hues show the constraint obtained 
in the 3 + 2 model and the 3 + 3 toy model introduced in section [TT| respectively. To obtain these 



results, we have used the methods described in section |III| to predict the ratio of the event rates 
at IceCube with and without sterile neutrinos, and we have then rescaled the published IceCube 
90% CL limit on the dark matter-nucleon scattering cross section crgg^STD (which was computed 
assuming standard oscillations) by this ratio. Specifically, if we denote the IceCube event rate by 
^STD, ^3+2 and A^3+3 for the standard oscillation, 3 + 2, and 3 + 3 scenarios, respectively, we 
compute the cross section limits in the 3 + 2 and 3 + 3 scenarios, 0-90,3+2 and 0-90,3+3, according to 

^3+2 

C^90,3+2 = cr90,STD^r7 , (15) 

YVsTD 

^3+3 /.^N 

C^90,3+3 = cr9o,STDT7 • [^^) 

YVsTD 

We see that the 3 + 2 model leads to a moderate weakening of the cross section limit, which 
can be understood from the fact that only electron neutrinos z/g and muon antineutrinos are 
substantially transformed into sterile states (see figure [T]), and that these transitions happen only 
for neutrinos with energies above several hundred GeV, whose contribution to the muon flux at 
IceCube is suppressed due to the large absorption probability in the Sun. In the 3 + 3 toy model, 
on the other hand, resonant flavor transitions happen already at lower energy (see figure [2]), and 
they happen for z/g, z/^ and Ur- 

As mentioned in section |ll| the effect could be even stronger if Amli, ^^51 5 and Amg^^ were 
negative (which might, however, require non-standard cosmology to be consistent). 



V. DEPENDENCE ON STERILE NEUTRINO PARAMETERS 



In section IV we have illustrated using two exemplary models how neutrino limits on dark 
matter capture and annihilation in the Sun are modified by oscillations into sterile neutrinos . We 
are now going to study more systematically how the worsening of these limits depends on the sterile 
neutrino parameters. We do this using the 3 + 3 toy model introduced in section |TT] since this model 
has only two new parameters (Og and Am^), but still covers the most important phenomenological 
aspects of more general sterile neutrino scenarios. 

We show in figure [6] the factor by which the IceCube limits on the spin-dependent dark matter- 
nucleon scattering cross section are weakened for a wide range of sin^ 29s and Am^ values. The 
shape of the contours can be understood as follows: At very large Am^, the new MSW resonances, 
equations ([9])-(ll) lie at a very high neutrino energy. For instance, at Am^ = 1 eV^, equation ^ 
yields a resonance energy of about 60 GeV at solar core densities, i.e. only neutrinos with E > 
60 GeV are affected by the resonance. Since very high energy neutrinos are mostly absorbed in 
the Sun, they do not contribute significantly to the IceCube limits. For somewhat lower Am^, 
the resonances move down in energy into the region relevant to IceCube. For too low Am^ or 
for too small ^5, on the other hand, MSW-enhanced flavor transitions become non-adiabatic (see 
equation ^ and related discussion), suppressing active-sterile transitions again. This happens 
first at high energy, which is why at low Am^ the correction factors shown in figure |6] are generally 
larger for dark matter annihilation into the soft bb channel than for annihilation in to the hard 
W^W- final state. 
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Dark Matter mass [GeV] 

Figure 5: IceCube limits on spin-dependent dark matter-proton scattering [27l [57] in scenarios with 
(black/gray dashed and dotted hnes) and without (black/gray sohd hues) sterile neutrinos compared to 
data from direct detection experiments [24l [63H65] and from Super-Kamiokande (SK) [25 (colored lines). 
Black lines correspond to results based on our Monte Carlo (MC) code, whereas gray lines are based on 
a simplified calculation which does not include secondary neutrinos (see text for details). We see that for 
the 3 + 2 scenario which provides the best fit to short baseline neutrino oscillation data, the limits are only 
moderately weakened. Our 3 + 3 toy model, on the other hand, illustrates that larger modifications are 
possible. 

VI. DISCUSSION AND CONCLUSIONS 

In this paper, we have shown how IceCube limits on dark matter capture and annihilation in 
the Sun are modified if eV-scale sterile neutrinos exist, as suggested by part of the short baseline 
oscillation data. Since IceCube is looking for high-energy neutrinos from dark matter annihilation 
in the center of the Sun, its results depend strongly on the oscillations of these neutrinos on 
their way out of the Sun. We have argued that in sterile neutrino scenarios new high-energy 
MSW resonances can lead to almost complete conversion of certain neutrino flavors into sterile 
states inside the Sun. In this case IceCube's constraints on dark matter-nucleon scattering can be 
significantly weakened, by a factor of two or more. 

This may have interesting implications if in the future dark matter is detected in a direct search 
or at the LHC, but the parameters determined there are in conflict with limits (or signals) from 
neutrino telescopes. If the allowed dark matter annihilation channels and branching fractions are 
established at the LHC, such a conflict could then provide a clear and strong hint for the existence 
of sterile neutrinos. With sufficient data, neutrino telescopes would even be able to contribute the 
determination of the active-sterile mixing parameters. 

Note added: While we were completing this work, reference [66j appeared on the arXiv, ad- 
dressing similar topics. 
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Flux ratio for 3 + 3 toy model, = 1 TeV 




sin2(26's) 



Figure 6: Weakening of IceCube limits on dark matter capture and annihilation in the Sun due to sterile 
neutrinos, assuming for illustrative purposes the 3 + 3 toy model introduced in section |TI[ The contours show 
the factor by which the IceCube limit on spin-dependent dark matter-proton scattering cross section for a 
1 TeV WIMP. Red dashed contours are for annihilation into (which yields a rather hard neutrino 

spectrum), blue solid contours are for annihilation into bb (which yields a much softer spectrum). At large 
Am^, oscillations into sterile neutrinos become less relevant because the active-sterile MSW resonances 
move to very high energies; at small Am^ or small sin^ 20s, the MSW transitions become non-adiabatic. 
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Appendix A: Numerics of neutrino oscillation probabilities 

In this appendix, we discuss the algorithm used to compute the neutrino oscillation probabilities 



in the Sun. As mentioned in section [III] we use an implementation of the well known Runge-Kutta 
(RK) algorithm [67j, namely the rkf45 algorithm implemented in the GNU Scientific Library [56j. 
In each iteration this algorithm uses a step function to evolve the neutrino state vector from a time 
to to a time to + At by approximately solving the Schrodinger equation, where At is chosen such 
that the optimal balance between speed and accuracy is achieved. Rather than working entirely 
in one basis, we transform the Schrodinger equation to an instantaneous interaction basis before 



13 



each step. This instantaneous interaction basis is defined by the transformation 

^Piit; to) = Sit, to) m = e^^°(*-*°V(t) , (Al) 
where the Hamiltonian has been separated in the fohowing manner 

H{t) = H{to) + AH{t;to), (A2) 

with 

Hoito) = ^UDU^ + Vito) , AHit; to) = Vit) - Vito) . (A3) 

Here V{t) is the neutrino matter potential (see equation Q), £^ the neutrino energy, U the leptonic 
mixing matrix, and D = diag(0, Am2i, Amg^^, ...). 

The Schrodinger equation in the interaction basis is 

= Hi^j (A4) 

with Hiit;to) = Sit,to) AH SHt,to). Since the matter potential changes slowly in the Sun and 
thus Hj is small, the RK algorithm can choose a larger step size At compared to a calculation in 
the flavor basis. Note that the elements Sjk of the transformation matrix S{t^ to) can be evaluated 
efficiently by computing Vjme~^^'^^^~^^\y^)mki where are the eigenvalues of Hq, and V is 
the matrix of the corresponding eigenvectors. After the evolution of the step has concluded we 
transform back to the flavor basis, 

i/j{to + At) = e-^^°(^°)^V/(^o + At) , (A5) 
and proceed to the next step, setting to ^ to + At. 
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